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Abstract. The 7r + 7r~ and qq annihilation mechanisms of dilepton pro- 

duction during relativistic nuclear collisions are studied. We focus on the mod- 
ifications caused by the specific features of in-medium pion and quark states 
rather than by medium modification of the p-meson spectral density. The 
main ingredient emerging in our approach is a form-factor of the multi-pion 
(multi-quark) system. Replacing the usual delta-function the form-factor plays 
the role of distribution which, in some sense, "connects" the total 4-momenta 
of the annihilating and outgoing particles. The difference between the c.m.s. 
velocities attributed to annihilating and outgoing particles is a particular con- 
sequence of this replacement and results in the appearance of a new factor in 
the formula for the lepton pair production rate. We obtained that the form- 
factor of the multi-pion (multi-quark) system causes broadening of the rate 
which is most pronounced for small invariant masses, in particular, we obtain a 
growth of the rate for the invariant masses below two masses of the annihilating 
particles. 
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1. Introduction 

The main objective in studying strongly interacting matter, which is formed during 
relativistic heavy ion collisions, is identification of the quark-gluon plasma. Lepton 
pairs and photons created during relativistic nucleus-nucleus collision do not interact 
with the highly excited nuclear matter, they leave the reaction zone without further 
rescattering. That is why, the dileptons (e + e~ and /i + pr pairs) observed in high- 
energy heavy ion collisions carry undistorted information on the dense early stage 
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of the reaction as well as on its dynamics. The enhancement with an invariant mass 
of 200 -j- 800 MeV observed by the CERES collaboration [Q] |2] in the production of 
dileptons has received a considerable attention (for the review, see Ref. [[3]). It was 
found that a large part of the observed enhancement is due to the medium effects 
(see Refs. [HJ [5] and references therein). Meanwhile, pion annihilation is the main 
source of dileptons which come from the hadron matter [[6l [7] (see also [H]). That 
is why, the proper analysis of the dilepton spectra obtained experimentally gives 
important data which probe the pion dynamics in the dense nuclear matter. The 
purpose of the present paper is to look once more on the 7r + 7r~ and qq annihilation 
mechanisms of dilepton production from the hadron and quark-gluon plasmas by 
accounting the medium-induced modifications of the dilepton spectrum. In order to 
do this, we concentrate on the modifications which are due rather to intramedium 
pion and quark states, than on the discussion of a modification of the p-meson 
spectral density. In accordance with our suggestions, the main features of a pion 
wave function follow from the fact that pions live a finite time in the system where 
they can take part in the annihilation reaction. As a consequence, the off-shell 
effects give an appreciable contribution to the features of the annihilation process 
specifically in the region of low invariant masses. Moreover, if pions are the entities 
of a local subsystem, then the spatial structure of the pion states is far from a 
plane-wave one and this also gives the essential contribution to the features of the 
dilepton spectrum. 

2. Annihilation of particles in finite space-time volume 

To carry out the outlined program, we assume that the pion liquid formed after the 
equilibration exists in a finite volume, and the confinement of pions to this volume 
is a direct consequence of the presence of the dense hadron environment which 
prevents the escape of pions during some mean lifetime r. The same can be assumed 
concerning a hot system of quarks which are confined to a quark-gluon droplet. 
So, we assume the system of pions (quarks) produced in high-energy heavy-ion 
collisions is effectively bounded in a finite volume. The pion (quark) wave functions 
cj)\(x), where A is a quantum number, satisfy the proper boundary conditions and 
belong to the complete set of functions. For instance, the stationary wave functions 
may be taken as the solutions of the Klein-Gordon equation (V 2 + k 2 ) 4>\ (x) = 
0, where k 2 — E 2 — m 2 , which satisfy the Dirichlet boundary condition on the 
surface S: (f)\(x)\s = 0. For the box boundary, we get </>k(x) = \/8/V Y\^ =1 0(Li — 
Xi)6(xi) sin (kiXi), where V = LiL 2 L 3 is the box volume, A = k = (fci, k 2 , fe), and 
components of the quasi- momentum run through the discrete set ki = irrii /Li with 
rii = 1, 2, 3, . . . For the spherical geometry, the normalized solutions are written as 
hlm(r) = e(R-r) (2/r) 1/2 J l+1/2 (kr)Y lm ($,a)/RJ l+3/2 (kR), where A = (k,l,m). 
Next, the field operators <f>{x) corresponding to the pion field should be expanded 
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in terms of these eigenfunctions, i.e. 

d 3 k 



(27r) 3 2w k 



[a(k)Mx)+b + Qt)<f>* k (x)] , 



(1) 



where a(k) and 6(k) are the annihilation operators of positive and negative pions, re- 
spectively. On the other hand, the states corresponding to confined particles can be 
written in a common way as </>k(x) 

III 



y/p(x)/V ( &k(x), where, for instance, yj p(x) 
9{R 



_ x [8 {Li — Xi)0{xi)] for a box and yj p(x) = 6{R — r) for a sphere, respectively. 
The function p(x) represents the information about the geometry of a reaction re- 
gion or cuts out the volume where the pions (quarks) can annihilate. Hence, for the 
evaluation of S- matrix elements wave functions (f>k{x) should be taken as the pion 
in-states once annihilating pions belong to finite system. The amplitude of the pion- 
pion annihilation to a lepton pair in the first non-vanishing approximation is calcu- 
lated via the chain 7r + 7r~ — ► p — ► 7* — * 11, where the p-meson appears as an inter- 
mediate state in accordance with the vector meson dominance. The matrix element 
of the reaction is (out|5|in) = — f d 4 xid 4 X2 (p+,p_ | T [HJ {x\) H\{x2)\ |ki,k2), 
where Hj{x) = —ejJ l {x)A^{x) and Hj{x) = —ejj i {x)A fJ, {x). It is remarkable that 
the pion density p(x) appears as a factor of the pion current. Indeed, 



3a(%) = -i<f{x) df, <p + {x) - 



V 



-i$(x) 8,, <l> + (x) 



(2) 



where the field operator $(a;) is defined in the same way as that in {T]) with the 
functions <frk{x) replaced by <&k(x). Because of this factorization, after the integra- 
tion over the vertex x the density p(x) automatically cuts out the volume, where 
the 7r + 7r~ annihilation reaction is running. At the same time, this means that the 
density p(x) determines the volume of quantum coherence, i.e. just the particles 
from this spatial domain are capable to annihilate one with another and make con- 
tribution to the amplitude of the reaction. To obtain the overall rate, it is necessary 
then to sum up the rates from every coherent domain of the fireball. 

For the sake of simplicity one can assume that the pion states can be approx- 
imately represented as 4>u{x) = \J p{x) /Ve~ tk ' x {p{x) is the 4-density of pions in 
the volume V where pions are in a local thermodynamic equilibrium). In essence, 
this approximation considers just one mode of the wave function ^(x) and reflects 
the qualitative features of the pion states in a real hadron plasma. In the frame of 
this approximation a simple calculation immediately shows that the S'-matrix ele- 
ment is proportional to the Fourier-transformed pion density p{x), i.e. (out|S , |in) cx 
p{k\ + hi — p + — p~), where k\ and ki are the 4-quasi-momenta of the initial pion 
states and p+ and p_ are the 4-momenta of the outgoing leptons. This means that 
the form-factor of the pion source p{k) stands here in place of the delta function 
which appears in the standard calculations, i.e. {2-k) a 8 a {K — P) — > p{K — P), 
where K = kx + ki and P = p + + p_ are the total (quasi-) momenta of pion and 
lepton pairs, respectively. An immediate consequence of this is a breaking down of 
the energy-momentum conservation in the s-channel of the reaction, which means 
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that the total momentum K of the pion pair is no longer equal exactly to the total 
momentum P of the lepton pair. The physical interpretation of this fact is rather 
obvious: the effect of the hadron environment on the pion subsystem which prevents 
the escape of pions from the fireball can be regarded during the time span r as the 
influence of an external nonstationary field. The latter, as known, breaks down the 
energy-momentum conservation. From now, the squared form-factor \p(K — P)\ 2 
of the pion system plays the role of a distribution which in some sense "connects" 
in s-channel the annihilating and outgoing particles instead of 5-function. Indeed, 
the number of produced lepton pairs from a finite pion system related to an 
element of the dilepton momentum space, reads 



(3) 



where 

5 4 (p++ P - -P)\A (k u k 2 ;p + ,p.)\ 2 . (5) 



d A p + d A p- 
(2n) 3 2E + {2irf2E^ 



Here, Ei = \Jm 2 + kf , i — 1, 2 for pions and E% = y/mf + pf , i = +, — for leptons, 
respectively. To obtain Eq.((3]), we represent the amplitude of the reaction as 

(outlS^lin) =p(k 1 + k 2 -p + -p-)Ao(k 1 ,k 2 ;p + ,p-). (6) 

Note that not only the form-factor p(K — P) contains information about the pion 
system. The amplitude carries also new important features, which are related to 
the violation of the energy-momentum conservation in the s-channel. The latter is a 
consequence of the medium effects through a partial confinement of the pion states 
inside fireball what results in the breaking of the translation invariance. Indeed, 
the pion-pion c.m.s. moves with the velocity vk = K/A'o, whereas the lepton- 
pair c.m.s. moves with the velocity vp = P/Pq. Hence, these two center-of- 
mass systems are " disconnected" now, that is why any quantity should be Lorentz- 
transformed when transferred from one c.m.s., to another. Reflection of this is 
the appearance of the correction factor 



ii m ( p - K r ! 



in the formula for the 



dilepton production rate (for details see [0 [TU]). 

By the broken brackets in Eq.([3]), we denote the thermal averaging over the 
pion quasi-momentum space. Actually, we assume a local thermal equilibrium 
in the multi-hadron (-pion) system. Hence, the Green's function, D^xi, 1E2) = 
($ + (x2)$(a;i)), which appears after thermal averaging, can be represented as 

D<( Xl ,x 2 )= J -0^e- lk - x A(k)f BE (k,X), (7) 
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where X = {x x + x 2 )/2, x = x x - x 2 , / BE (fc, X) = {exp [/3(X)(k ■ u(X) - n(X))] - 
is the Bose-Einstein distribution function, which depends on space-time vari- 
ables, X = (X°,X), (3 is the inverse temperature, u(X) is the hydrodynamical 
velocity and p is the chemical potential. For the ideal gas (infinite life time of the 
system) the spectral function A(k°,'k) indicates that all states are on the mass- 
shell: A(k) = 4ir5(k 2 — m 2 )9(k°). In the interacting system the spectral function 
reflects a collision broadening of the states which includes as well a global decay of 
the system if collisions in the system exist during finite time span. For instance, 
the fireball, which is nothing else as a system of strongly interacting particles, lives 
until its decay, i.e. starting from the creation till the freeze-out, after which there 
are no strong interactions between particles at all. In our further consideration we 
will take into account just a global decay of the multi-pion (multi-quark) system. 
Assuming a proper model of the spectral function, A(k), which responsible in the 
present approach for finite life time of the system, we incorporate it together with 
the spatial density p(x) to the global system form-factor p(x). 

Concerning the physical meaning of Eq. |3]), we note that one can regard it as 
the averaging of the random quantity ( d A< Kd i p ) w ith the help of the distribution 
function \p(K — P)\ 2 centered around the value P, which is fixed by experimental 
measurement. In this sense, the hadron medium holding pions in a local spatial 
region for some time, which is expressed as the local pion distribution p{x), plays 
the role of an environment randomizing the pion source. This randomization is 
a purely quantum one in contrast to the thermal randomization of the multi-pion 
system which is already included to the quantity ( grf^Fp )- 



3. Dilepton emission rates 



In order to transform the distribution of the number of created lepton pairs over 
the dilepton momentum space to the distribution over invariant masses, one has 



to perform additional integration using {dN^ /d 4 P) from i.e. 



J 2Pn 



dFP 



where Pq = \J M 2 + P 2 . This results in: 



dM' 2 
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2m 2 
~M 2 



\F„(M 2 )f 



1 



[P-K) 
M 2 K 2 



1 



d 3 P 
2Po" 



d 4 K 



M 2 



(8) 



where we take the Boltzmann distribution Jbe(E) ~ exp(— j3E). Note, that during 
integration with respect to a 4-momentum K one should keep the invariant mass 
of a pion pair, = \J K 2 , not less than two pion masses. On the other hand, 
possible finite values of the distribution (^^jp-^ below the two-pion mass threshold 
can occur just due to the presence of the pion system form- factor p(K — P). The 
factor in the square brackets on the r.h.s. of ijHJ) is a correction which is due to the 
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Lorentz transformation of the quantity (ki — k2) 2 from the dilepton c.m.s. to the 
pion-pion c.m.s. This factor gives a remarkable contribution to the dilepton spec- 
trum for invariant masses below the two-pion mass value. Its influence is especially 
pronounced for e + e~ production as was shown in [[9l I10j . 

In Eq.® the p-meson form-factor, F 7T (M 2 ), is a vacuum one. Actually, there are 
two ways to take into account effects of the hadron medium: first, one can account 
for distortion of the pion states caused by dense environment; second, one can look 
for p- meson polarization effects during its passing through the hadron environment. 
In the present paper we choose the first way of accounting for the medium effects (see 
also [HI]). Just to elucidate as much as possible the consequences of the contraction 
of the pion states in the hadron medium and to prevent a double counting we take 
the vacuum p-meson form-factor. 

For particular evaluations we take as a model of the pion system the Gaussian 
distribution of the particles in space and the Gaussian decay of the system (this 
form- factor succeeded in HBT intcrfcrometry) : 




p(x) — exp 



x 2 - 2(u ■ xf 
2R 2 



(9) 



where u is the hydrodynamical (collective) veloc- 
ity of the element of the total system which is in 
a local thermodynamic equilibrium; R is the spa- 
tial radius of the element. In the rest frame of 
the element the form- factor looks like: po(t,v) = 
exp [— (t 2 + r 2 )/2i? 2 ]. To get a proper interpretation 
in terms of the mean life time of the system element, r, 
one needs to make a scale transformation during inte- 
gration over the time variable: dt po (t, r) F(t, r) = 
dt exp (-t 2 /2r 2 - v 2 /2R 2 )Fh{t, r), where F (t, r) 
= ZF(ft,r). 

Meanwhile, it can be another choice of the pion 
source function. Indeed, one can choose, for instance, 
a geometry with sharp boundaries which are determined by the ^-function. To 
show that the final answer is not sensitive to the form of the cutting function we 
compare two form-factors (normalized to the unit volume) which correspond to 
the Gaussian distribution p(r) — exp (— r 2 /2R 2 ) and to the ^-function distribution 
p(r) = 6(R~ |r|) (see Fig. 1). Only a slight difference between these form-factors is 
seen and, therefore, the choice of pion source distribution does not affect much the 
dilepton production rate. 



Fig. 1. Comparison of the 
Gaussian and 8— function form- 
factors |p(q)| 2 /V. 



We evaluate the rate in the rapidity window, y n 
sponds to CERES experimental conditions [[TJ [2] 



< v < 



y s yn 



which corre- 
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dN 
d 4 xd 4 P 
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where tanhy = P 3 /P°, 



duction rates dR^+ e _ /dMdy and dR (l + _/dMdy for electron-positron and muon- 
muon pairs, respectively, in pion-pion annihilation arc depicted in Fig. 2. Note, 
the calculations are carried out in the frame of the element of the system where 
particles are in a local thermal equilibrium. Different curves correspond to the dif- 
ferent "spatial sizes" R and different "lifetimes" t (for particular values of these 
parameters see Fig. 2) of a hot pion system at the temperature T = 180 MeV. 
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(P 1 ) 2 + (P 2 ) 2 . The results of evaluation of the pro- 
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Fig. 2. Rates of electron-positron (left panel) and muon-muon (right panel) productions in pion- 
pion annihilation in a small finite system, T = 180 MeV. 

For comparison, we present in Fig. 3 the results of evaluation of the rate 
dR^ p } _ /dMdy of electron-positron pair production in quark- ant iquark annihilation 
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Fig. 3. Electron-positron production rates in quark- ant iquark annihilation in a hot QGP drop. 

in a hot QGP drop, T — 180 MeV. The evaluation was carried out in the frame of 
the quark drop under the same assumptions as for pion-pion annihilation. As in the 
previous case, an increase in the rate with decrease in the invariant mass up to two 
electron masses is seen. This real threshold is close to the total mass of annihilating 
quarks M = 2m q % 10 MeV/c 2 . 
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4. Conclusions 

We notice that the production rate in a finite small system differs from the rate in 
an infinite pion gas (solid curve) where pion in-states can be taken as plane waves. 
The deviation bigger when the parameters R and r are smaller. Of course, this 
is a reflection of the uncertainty principle which is realized by the presence of the 
distribution \p(K — P)\ 2 as the integrand factor in ([3|). Basically, the presence of 
the form-factor of the multi-pion system will result in a broadening of the rate for 
small invariant masses M ^ 800 MeV/c 2 which is wider at the smaller parameters 
R and r. This seems natural because the quantum fluctuations of the momentum 
are more pronounced in smaller systems. We emphasize as well that the behavior of 
the curves in Fig. 2 which correspond to a finite system has a similar to the CERES 
data tendency [ HI E] • 

The same behavior of the rate is seen in a hot quark drop (see Fig. 3): for small 
parameters R and r in the region of small invariant masses M ^ 500 MeV/c 2 , as 
compared to the rate for infinite parameters R = oo, r = oo, is due to a rise of 
quantum fluctuations which are evidently bigger for a smaller size of the QGP drop. 

Note that the enhancement of the dilepton production rate for the low invariant 
mass region is much more sensitive to the variation in the spatial size of a many- 
particle (pion, quark) system than to that in the system lifetime (see Fig. 3, right 
panel) [[T2]. 
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